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ABSTRACT 

We examine the effects of turbulent intermittency on the deflagration to det- 
onation transition (DDT) in Type la supernovae. The Zel'dovich mechanism 
for DDT requires the formation of a nearly isothermal region of mixed ash and 
fuel that is larger than a critical size. We primarily consider the hypothesis by 
Khokhlov et al. and Niemeyer and Woosley that the nearly isothermal, mixed 
region is produced when the flame makes the transition to the distributed regime. 
We use two models for the distribution of the turbulent velocity fluctuations to 
estimate the probability as a function of the density in the exploding white dwarf 
that a given region of critical size is in the distributed regime due to strong local 
turbulent stretching of the flame structure. We also estimate lower limits on the 
number of such regions as a function of density. We find that the distributed 
regime, and hence perhaps DDT, occurs in a local region of critical size at a 
density at least a factor of 2 — 3 larger than predicted for mean conditions that 
neglect intermittency. This factor brings the transition density to be much larger 
than the empirical value from observations in most situations. We also consider 
the intermittency effect on the more stringent conditions for DDT by Lisewski et 
al. and Woosley. We find that a turbulent velocity of 10^ cm/s in a region of size 
10^ cm, required by Lisewski et al., is rare. We expect that intermittency gives 
a weaker effect on the Woosley model with stronger criterion. The predicted 
transition density from this criterion remains below 10'' g/cm^ after accounting 
for intermittency using our intermittency models. 

Subject headings: stars: interiors — supernovae: general — turbulence 



Introduction 



A successful model for Type la Supernova (SNIa) explosions is required to produce a de- 
flagration to detonation transition (DDT) by observational constraints. A pure deflagration 



- 2 - 



model gives exploding kinetic energy lower than observed (Khokhlov 1991; Gamezo et al. 
2003; Ropke & Hillebrandt 2005) and pure detonation leads to overproduction of iron group 
elements and too little intermediate elements (Branch et al. 1982, 1983). The density ptr 
at which the transition occurs determines the amount of the nickel produced (Hofiich 1995; 
Hoflich and Khokhlov 1996; Dominguez, Hofiich & Stranicro 2001). Therefore a prediction 
of ptr-i consistent with the observed nickel production, is essential to a DDT theory for SNe 
la . 

The mechanism by which the DDT occurs still remains a mystery The most studied 
candidate is the Zcl'dovich mechanism, which requires the existence of an almost isother- 
mal region of mixed ash and fuel that is larger than a critical size to drive a supersonic 
shock that is sufficiently strong to sweep over the entire star (Khokhlov et al. 1997, here- 
after KOW; Niemeyer and Woosley 1997, hereafter NW). One hypothesis is that a nearly 
isothermal region is produced by turbulent preconditioning. KOW argued that, to produce 
an almost isothermal mixture of ash and fuel, the laminar flame must be quenched by tur- 
bulent stretching, at least locally. This might allow the cold fuel to mix with the ash both 
thermally by electron conduction and chemically by diffusivity without being burned. They 
assumed that the criterion to quench a flame is that the turbulent velocity at the laminar 
flame thickness must be larger than the laminar flame speed. NW gave a similar argument 
based on the distributed flame burning regime in turbulent combustion. The criterion for a 
distributed flame is expressed in terms of the Gibson scale at which the turbulence velocity 
equals the laminar flame speed. If the Gibson scale is smaller than the laminar flame thick- 
ness, turbulent stretching can generate structures within the flame and the flame is in the 
distributed regime. NW speculated that in this regime flames can be temporally quenched 
in some regions, which can host the detonation after being homogenized in temperature and 
composition by turbulent mixing. The criterion for the distributed regime is equivalent to 
that for flame quenching used by KOW. Both criteria give the same condition on the turbu- 
lence intensity for given laminar flame properties (see §2). As the density in the star drops 
due to the overall expansion, it is easier for turbulence to affect the laminar flame because 
of the decrease in the flame speed and the increase in flame thickness. With presumed tur- 
bulence parameters, the criterion for the turbulence intensity, determined by the robustness 
of laminar flames disturbed by turbulent motions, translates to a transition density ptr for 
the DDT. 

Several uncertainties exist in the simple model given by these two early studies. First, 
it is not clear whether the criterion used by KOW, equivalent to that for a distributed 
regime (NW), is sufficient for flame breaking. How, or even if, flames are quenched is still an 
open question. Second, it is uncertain whether (local) flame quenching is indeed necessary 
to produce a nearly isothermal region. Finally, later studies by Lisewski, Hillebrandt and 
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Woosley (2000) (hereafter Lisewski et al 2000; see also Lisewski et al. 2000b) and Woosley 
(2007) find that entering the distributed regime, while probably a necessary condition, is 
not sufficient for the DDT to occur. Based on a requirement for turbulent transport to be 
efficient at producing a shallow temperature and composition gradient around the laminar 
flame, Lisewski et al. (2000) find that the turbulent velocity at the scale 10^ cm needed for 
a detonation is very large, ~ 10^ cm/s. Woosley (2007) claims that the DDT occurs only 
when the turbulent flame thickness exceeds a critical length scale. We show in §2 that the 
two criteria, although arising from different physical considerations, are basically equivalent. 
The corresponding condition is more stringent than that assumed by KOW and NW. 

In this paper we examine the effect of turbulent intcrmittcncy on the onset of distributed 
burning that may relate to the DDT. Despite the uncertainties listed above, we will mainly 
consider the model by KOW and NW and use it to illustrate the potential importance of 
intermittency in SN la explosions. Our calculations can be applied to the criteria by Lisewski 
et al. (2000) and Woosley (2007) in a straightforward way. A quantitative analysis using 
their criteria requires data for laminar flame properties and critical length scales at densities 
below 10'' g cm~^ that are not immediately available (see §2). We give a qualitative discussion 
of the intermittency effect on their DDT models. 

Intermittency is an important concept in turbulence theory. It is characterized by intense 
local events, e.g., strong stretching at small scales, which occur at a frequency much larger 
than predicted from a Gaussian distribution (see, e.g., Frisch 1995). The physical origin of 
intermittency in turbulent flows is the spatial inhomogcncity in the energy dissipation rate: 
most kinetic energy is viscously dissipated in the finest structures, e.g., vortex tubes, which 
occupy only a small volume fraction. These rare but intense dissipative structures give rise 
to a spatially inhomogeneous and intermittent distribution for the turbulent intensity and 
the stretching rate. Intermittency is shown as broad exponential tails in the probability 
distribution for the stretching rate or the dissipation rate at small scales (see §3). The tails 
get broader at smaller scales, meaning that the probability of finding an extreme turbulent 
stretching rate or intensity increases with decreasing scales. 

According to ID simulation results by KOW and NW, the critical size, Ic, of the isother- 
mal region required for a DDT via the Zcl'dovich mechanism is much smaller, especially at 
large densities, than the expected integral length scale for the buoyancy-driven turbulence 
in SN la explosions. This suggests that only a small flame region with a sufficiently strong 
local turbulence intensity may be needed to trigger a detonation. Turbulent intermittency, 
which indicates the existence of regions of small sizes where the turbulent stretching is much 
larger than the average value over the flow, is therefore expected to have important conse- 
quences for DDT. The transition could happen earlier at a higher transition density ptr than 
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predicted by models using the average turbulent intensity. At higher densities, much larger 
turbulent intensity is required for the DDT, but the rapid decrease of Ic with increasing den- 
sity makes an earlier DDT possible for two reasons. First, the probability is larger to find 
regions of smaller sizes Ic with extreme turbulent stretching rate or intensity. Second, there 
are more regions of smaller size available as candidates to host the detonation. Clearly, the 
intermittency effect accounts for the intuitive dependence of ptr on 1^. the smaller the critical 
size, the easier it may be for the transition to happen. To what degree the intermittency 
effect increases ptr is the main question we investigate in this paper. 

In §2, we review the criteria for the DDT in models by KOW, NW, Lisewski et al. 
(2007) and Woosley (2007) and formulate a new criterion taking into account the effect of 
intermittency. We describe two intermittency models by Oboukhov (1962) and Kolmogorov 
(1962) and by She and Leveque (1994) in §3. Using the intermittency models, we evaluate 
the transition density from the new criteria in §4. Our results are summarized and discussed 
in 55. 



2. Criteria for the DDT 

The criterion used in NW forjudging whether a flame is in the distributed regime, which 
was also assumed to be the condition for the DDT, is to compare the Gibson scale Iq with the 
laminar flame thickness //. The Gibson scale is defined such that Su{Ig) = Si where 6u{l) is 
the amplitude of the velocity fluctuations at the scale / (or equivalently the velocity difference 
over a scale /, i.e., 6u{l) = u{l + x) — u{x)) and Si is the laminar flame speed 0. If h ;^ If, 
the turbulence cannot internally disturb the flame and the turbulence in effect wrinkles the 
flame. This is called the flamelet regime. Only when Ig ^ //, can turbulence stretch the 
flame efficiently to generate structures within the flame and the turbulent combustion enters 
the distributed regime. The condition Iq ^ If is equivalent to 5u{lf) ^ 5u{Ig) = Si since 
5u(l) is an increasing function of the scale H The latter, which means that the turbulent 



^Note that NW, accounting for the cellular stabilization effect against instabilities, e.g., the Landau- 
Darrieus instability, defined Iq as the scale where the turbulent velocity exceeds the effective cellular flame 
speed. This does not introduce a significant difference in the estimate of la since the effective cellular speed 
is close to the laminar speed and has a very weak dependence on scale; see their Fig. 1. 

^This is also equivalent to the diffusivity criterion by Niemeyer and Kerstein (f997) for the onset of 
distributed regime and flame extinction at Prandtl number larger than unity, which is the case for white 
dwarfs. Their criterion was motivated by the observation that, at Prandtl number different from unity, two 
previous criteria proposed for the flamelet breakdown and for the flame quenching, using the ratio of the flow 
viscous length scale to the flame thickness and the ratio of the viscous timescale to the reaction timescale, 
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velocity fluctuation 5u{lf) at the scale of the flame thickness If is larger than the laminar 
flame speed, is the criterion used in KOW for flame quenching and the DDT. 

Following KOW, we introduce a factor of X ~ 1 in the criterion to account for the 
uncertainty in the flame breaking mechanism, i.e., Su{lf) > KSi. We will consider two 
values for K, i.e., K — 1 and K — % (KOW; For K — to quench a flame, the Gibson 
scale has to be = 512 times smaller than the flame width). This criterion can also be 
written in terms of timescales. Noting that the turbulent stretching timescale, r^, at the 
flame thickness is rt{lf) = lf/6u{lf) and that the nuclear reaction timescale, r„, is related 
to the flame speed r„ = If/ Si, the criterion is equivalent to rt{lf) < Tn/K, i.e., to break 
the flame the stretching timescale at the flame thickness must be smaller than the nuclear 
burning timescale (see Niemeyer and Kerstein 1997). 

To apply this criterion, the Kolmogorov (1941) scaling 5u(l) = e^/^/^/^ is usually used to 
calculate 5u{lf ) from the turbulent velocity fluctuations at large scales where e is the average 
dissipation rate in the flow. From this scahng, the criterion can be written as (KOW, NW), 

eV3/y3 ^ 

or 

e > K'Sf/lf = Khf (2) 

where cj is defined as Sf /If. Although we use the convenient criterion (2) in terms of the 
dissipation rate in our calculations, the turbulent stretching is more fundamental and we 
will use the concept of the flame stretching in our discussions. 

The laminar flame speed and thickness depend on the chemical composition and the 
density (Timmcs and Wooslcy 1992, KOW). In Table 1, we list the fiamc speed, the flame 
thickness as a function of density for a white dwarf with half carbon and half oxygen, mainly 
taken from Timmes and Woosley (1992). The laminar speed decreases and the thickness 
increases quickly with decreasing density p, therefore e/ decreases rapidly with decreasing 
p as shown in Table 1. The average dissipation rate is estimated to be e = U^/L where 
U and L are the characteristic velocity and length scales of the turbulence, normally set 
by motions on the large, driving scale. At large scales, the turbulence is driven by the 
Rayleigh- Taylor instability. The length scale L might be expected to be about the size, Rf, 
of the fiame region, L Rf 10^ cm and the velocity scale to be about the Rayleigh- Taylor 
velocity at this scale U ~ ^.h^fg^fjL ~ 10® cm/s where the effective gravity is taken to be 
~ g^ff = 5 X 10® cm/s^ (KOW, NW). Khokhlov (1995), however, showed that motions at 



respectively, are not equivalent. 
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scales larger than 10^ cm freeze out due to the overall expansion of the star. In that case, 
L ~ 10^ — 10^ cm and U ~ 10^ cm/s. We will take U and L as parameters. Note that the 
criterion eq (2) depends on U and L through the dissipation rate. Given the dissipation rate 
e, the critical density below which the inequality (cq 2) is satisfied can be obtained using 
€/ as a function of p in Table 1. For example, if [/ ~ 100 km/s and L ~ 100 km, e ~ 10^^ 
cm^/s^ and we find that, from interpolation in Table 1, e is larger than ej at a density less 
than ~ 4 X 10''' g/cm^. Therefore criterion (2) predicts a transition density p^^ — 4 x 10^ 
g/cm^ lor K — 1 (see KOW and NW). II K — the predicted transition density is smaller, 
Ptr ~ 1.5 X 10^ g/cm^. In the second hne of Table 2, we give the predicted ptr for different 
values of the parameters, which decreases with decreasing dissipation rate e. The numbers 
in parenthesis correspond io K — 9). 

When using the criterion eq (2), we need to keep in mind that the spatial fluctuations 
of e (see §3) are completely neglected and the criterion only applies to the overall situation 
in the combustion flow. We will refer to this criterion as the mean criterion. When the mean 
criterion is met, the only implication is that the combustion is in the distributed regime 
in general. Considering the intcrmittcncy of turbulence, i.e., the spatially inhomogeneous 
distribution of the stretching strength, there can be places where the stretching rate is 
much weaker than the average. These places could still be in the fiamelet stage while most 
other places are in the distributed regime. Or conversely, even if the mean criterion (2) is 
not satisfied, one cannot exclude the possibility of there existing a region that experiences 
strong stretching and gets into the distributed regime when most of the structure is still in the 
fiamelet regime. This latter fact is important for the deflagration to detonation transition. 
The fact that the DDT does not require the entire star to be in the distributed regime but 
instead only needs a region of size much smaller than the white dwarf radius (see below, 
KOW), coupled with the intrinsic intcrmittcncy, suggests that DDT could occur earlier than 
predicted by cq (2) and hence at a larger transition density. The detonation can be triggered 
locally when a region appears that is larger than the critical size and enters the distributed 
regime due to a strong local stretching. It is important to study the degree to which this 
intermittency effect increases the transition density, which is constrained by observations. 
Clearly the answer depends on the critical size, which we consider next. 

The question of how large the isothermal region with well-mixed ash and fuel has to be 
for a detonation was studied by KOW (see also NW). In their model, the DDT occurs via the 
Zel'dovich mechanism (Zel'dovich et al. 1970) where the mixed region begins spontaneous 
ignition at the place with the minimum induction time, and the flame propagates with a 
phase speed equal to the inverse of the spatial gradient of the induction time, which is large 
for nearly isothermal and well-mixed regions and is not limited by the speed of sound. As 
the phase speed decreases below the Chapman- Jouget speed, a shock forms just ahead of the 
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flame front. Whether this shock can explode the whole star depends on the strength of the 
shock when entering the pure fuel, which is determined by the size of the isothermal region. 
If the isothermal region is small and the shock is weak, the flame front and shock separate 
with the flame front lagging behind the shock and the shock cannot make the whole star 
explode. The critical strength of the shock corresponds to a critical size of the isothermal 
region, over which the shock can be strengthened. Using ID simulations, KOW and NW 
obtained the critical size, Ic, which depends on the density and the chemical composition. 
It is interesting to note that at early time when the density is large, the required size is 
much smaller than that at later times. We will show this has important consequences. The 
critical size is much smaller than what current numerical simulations can resolve, therefore 
the problem of the intermittent stretching at scale Ic cannot be addressed by simulations. 

We need a local criterion to check whether a region of a given size I, in particular Ic, is 
in the distributed regime or not. For that purpose, we use a local average dissipation rate ei 
(see eq. 8 in §3 for a definition) in a region of size / to replace e. Following the same argument 
that leads to equation (2), the criterion for a region of size I being in the distributed regime 
is 

6, > K'ef (3) 

where we have used the refined similarity hypothesis by Kolmogorov (1962) (see eq. 9 in 
§3). Due to the random nature of turbulent fiows, e; is stochastic, and a statistical approach 
is necessary. We therefore ask the question: what is the probability that any region of size 
/ is in the distributed regime? This is given by the cumulative probability P(e; > K^ef). 
To answer this question, we need the probability distribution P(e/) of e^. Fortunately, this 
distribution has been extensively studied in the intermittency models for turbulence, which 
we describe in §3. Although these models were originally proposed for homogeneous and 
isotropic turbulence, we will assume they apply to SNe la where the turbulence is stratified 
and may not arrive at isotropy even at very small scales. Once the distribution is specified, 
one can calculate the probability of finding that a region of given size Ic is in the distributed 
regime, 

oo 

P{ei^>Khf)^ J P{eiMei. (4) 

K^ef 

which depends on the density through Ic and e/. An immediate examination of eq (4) shows 
that, at larger density, the lower limit of the integral K^ef is larger because of the fast flame 
speed and the small flame thickness. This tends to decrease the probability. However, at 
larger density, is smaller and the intermittency of turbulence tells us that the tail of the 
distribution P(e/) is broader for smaller /. This tends to counteract the decrease of the 
cumulative probability due to the larger lower integral limit at higher densities. 
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Purthermore, for smaller Ic, there are more regions of size Ic available in the star. This 
could make the transition occur significantly earlier with a transition density considerably 
larger than predicted by eq (2). We need to multiply the probability that a given region of 
size Ic is in the distributed regime by the number, A^^^^, of regions of size Ic available in order 
to calculate the number of regions that are both larger than Ic and in the distributed regime 
at any given density. We assume that the deflagration to detonation transition happens 
when 

Ni^ X P{ei^ > Khf) = 1. (5) 

Since we are concerned with the flame being stretched into the distributed regime, only 
locations around the flame front are of interest when calculating Ni^. Therefore, we only 
count regions in the vicinity of the flame front. Ni^ depends on the size, Rf,oi the flame 
region and the flame geometry. A typical value for Rf is 10^ cm (Khokhlov 1995), which 
could be smaller at an earher time. We will set i?/ ~ L in our calculations in order to decrease 
the number of parameters. Note that Rf > L when the freezeout effect is considered and 
therefore the number A*";^ we use is a lower limit. If the flame region is a 2D spherical front, 
A^";^ ~ AnR'j/ll. If the flame structure is highly convoluted, it may have a fractal dimension 
larger than 2. In that case, Ni^^ is larger. The upper limit for A^;^ is ^ 47ri?j/3Zj?, which 
applies if the flame geometry is close to 3D. Again, we take the lower limit A^;^ = AnR'j/ll, 
thus the transition density we will get is a lower limit. 

As discussed in the Introduction, Lisewski et al. (2000) and Woosley (2007) find that 
entering the distributed regime is not sufficient for the DDT to occur and give criteria 
stronger than that used in KOW and NW. Lisewski et al. (2000) considered how turbulent 
transport affects the temperature and composition profile around a laminar flame. They 
assumed that, at any point, turbulence translates the temperature and composition by a 
distance It, over which turbulence can transport during a local induction time Tj. The 
distance If is a function of position since depends on local temperature and composition. 
It is estimated by the length scale of a turbulent eddy with turnover time equal to Tj, 
i.e., lt/5u(lt) = Ti. Using the Kolmogorov (1941) scaling, we get k = e^^^rf^"^. For given 
turbulence intensity, temperature and composition profiles around a laminar fiame front can 
be calculated from the translation. Clearly, more efficient turbulent transport gives shallower 
temperature and composition profile, which is needed for detonation. By checking whether 
the resulting profiles, as initial conditions to solve the ID hydrodynamic equations, can lead 
to a detonation, Lisewski et al. (2000), obtained a condition for the DDT on the turbulent 
intensity. They found that, for a successful detonation, the turbulent velocity has to be 
~ 10^ cm/s at the scale 10^ cm. This condition is stronger than just entering the distributed 
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regimqj. Since the expected turbulent velocity at scale 10^ cm is 10^ — 10^ cm/s, Lisewski 
et al. (2000) concluded that a DDT via the Zel'dovich mechanism in SNe la is unlikely. 
However, considering the spatial inhomogeneity of turbulent intensity, i.e., intermittency, it 
is possible for regions of size 10^ cm with large enough turbulent velocity to arise. 

The result of Lisewski et al. (2000) motivated Ropke (2007) to study the probability of 
finding a region of size 10^ cm with a turbulent rms velocity of ~ 10^ cm/s. Using data from 
3D numerical simulations with a turbulent subgrid-scale method, Ropke (2007) analyzed the 
velocity fluctuations at the grid size (10^ cm) and obtained a fat exponential tail for large 
velocity fluctuations that extends up to 10^ cm/s. The large velocity fluctuations seem likely 
to be located at the trailing edge of a bubble-like feature (Ropke 2007). This confirms the 
intermittency in the turbulent combustion flow in SNe la; there exist grid cells where the 
turbulent intensity is much stronger than the average. From the probability of finding a 
grid cell with required turbulent intensity, Ropke concluded that the DDT triggered by a 
local cell with large velocity fluctuations is possible but probably rare. In our notations, 
the probability is given by -P(eio6cm > 10^*^ cm^/s^) where 10^^ cm^/s^ corresponds to the 
dissipation rate in a region of size 10^ cm with a rms velocity of 10® cm/s. We will calculate 
this probability and consider the availability of such regions using two intermittency models 
given in §3 and compare with the results of Ropke (2007) in §4. 

Woosley (2007) proposed a new criterion for the DDT based on a calculation of the 
distributed flame width using an eddy diffusivity approximation. Making an analogy to 
the estimate of the laminar flame thickness, Woosley (2007) obtained the distributed flame 
width A from the equation A ~ (i5(A)r„)^/^ where -D(A) = Su{X)X is the eddy diffusivity 
at scale A and r„ is the nuclear reaction timescale. Using the Kolmogorov (1941) scaling 
for 6u{X), the distributed flame width is given by A = e^/^Tn''^ (note that this formula for 
A is similar to It in Lisewski et al. (2000)). Woosley (2007) assumed that the condition for 
detonation is that the minimum burning timescale in the distributed flame is smaller than 
the sound crossing time over the distributed flame width A, or equivalently, A ^ r^^nic where 
"^sonic is sound crossing length over the minimum burning timescale in the distributed 
flame. The minimum sound crossing length is thus the critical size of the distributed flame 
width for detonation. The criterion A ^ r^^nic equivalent to e {''^^^niJ^/'^n- Noting that 



■^This condition can be converted into a form that can be directly compared with eq. 2. Roughly speaking, 
the physical condition for a detonation in this model is that It at the laminar flame front is larger than Ic, 
i.e., a shallow temperature gradient can be produced over a critical size around the flame front. Requiring 
k > Ic, we get the criterion e > {lc/lf)'^{Tn/Ti)^ef. Considering that n defined by Lisewski et al. (2000) is 
smaller than the nuclear timescale t„ and that Ic 3> If, this condition is much stronger than the condition 
of eq. 2. Note that this condition is similar to the criterion of Woosley (2007) given below. 
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Tn — If /Si, the condition can be written as e ^ {''^TonidhY'^f- Since r^^^ given in Table 4 
of Woosley (2007) is close to Ic listed in Table 1 in this paper, we will use Ic instead of r^^^ 
for simplicity, i.e., 

e ^ ilc/lffef. (6) 

which is much stronger than eq (2) because 1^ is much larger than the laminar flame thickness 
If. This condition can be used to determine the transition density Pt,- by the same calculation 
process as in the case of the criterion cq (2). Note that, except a factor of (T„/rj)"^, this 
criterion is basically equivalent to that given in footnote (3) for the requirement by Lisewski 
et al. (2000). We find that, for the reasonable turbulence parameters listed in Table 2, the 
criterion results in a transition density below 10^ g cm~^ and we cannot specify it due to the 
lack of data for 1^ Si and If at densities below 10^ g cm~^. 

In his estimate for ptj., Woosley (2007) used U = 10^ cm/s at scale L — 10^ cm through- 
out the calculations, based on the result by Ropke (2007) on the possibility of the existence of 
regions of size 10^ cm with a rms velocity of 10^ cm/s. With these turbulence parameters he 
derived ptr — 10^ g cm~^. Clearly, in Woosley 's calculation, the intermittency effect implic- 
itly contributes to the transition density obtained because, as discussed earlier, a turbulent 
rms velocity of 10® cm/s at 10^ cm can only arise from intermittency. 

The intermittency effect for the criterion of Woosley (2007) can be included more con- 
sistently in our formulation. Instead of considering a single special scale 10^ cm, our model 

specifics intermittency over a continuous range of scales corresponding to critical sizes at 
different densities. Following the same steps that lead to eq (5), we incorporate the inter- 
mittency effect in the DDT model of Woosley (2007) and obtain a criterion, 

Ni^ X P(q, > {lo/lffef) = 1. (7) 

which only differs from eq. 5 by the lower limit in the cumulative probability. We will discuss 
about this criterion in §4. 

We point out that the eddy diffusivity method used by Woosley (2007) to approximate 
the combined action of the turbulent advection and the microscopic diffusivity is an oversim- 
plification. This procedure implicitly assumes a smooth structure in the distributed fiame 
and neglects the fiuctuations of temperature and concentration, which may be important in 
determining the effective width of distributed flames. 

3. Intermittency 

Kolmogorov's 1941 theory assumes that the energy transfer in the inertial range is equal 
to the average dissipation rate e in the flow and is the same throughout the inertial scales 
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down to the viscous scale where the kinetic energy is removed. This assumption, together 
with the similarity hypothesis, predicts that the statistics of the velocity difference (or the 
velocity fluctuations) at any inertial scale is completely determined by the average dissipation 
rate e. However, fluctuations in the dissipation rate clearly exist as can be seen from the 
formula for the local viscous dissipation rate, e(x, t) = | JZi^i'^^j + djUiY, which is a function 

of the fluctuating velocity field. The spatial fluctuations in e are well-illustrated by the intense 
dissipation structures at small scales such as vortex tubes. This effect needs to be taken into 
account for a more accurate prediction of the scaling behavior of the velocity difference 
(Landau and Lifshitz 1944). The statistics of the velocity difference over a separation / 
depends on the distribution of the dissipation rate over regions of size I, which is defined as 
(e.g., Kolmogorov 1962), 

ei{^,t)^^ j e(x + x',i)dx'. (8) 

|x'|</ 

Clearly, the mean of e; is equal to e and thus is independent of This means that the 
average energy flux over all the inertial scales is constant. The e; distribution is essential to 
the intermittency models for turbulence. Note that this distribution is exactly what we need 
in our calculations for the transition of the turbulent combustion to the distributed regime 
by turbulent stretching and quenching described in §2, eqs (4) and (5). 

Intermittency in turbulence is usually expressed in terms of the scaling behavior of the 
structure functions {5u{iy) ~ where 5u{l) — u{x + 1) — u{x) is the (longitudinal) velocity 
difference and Cp is the scaling exponent for the pth- order structure function. Kolmogorov's 
1941 theory predicts that the exponent goes with p as = p/3. However experimental 
data (e.g., Anselmet et al. 1984) have shown departure from this linear relation and (p 
increases significantly slower than p/3 at large p. This "anomalous" scaling is referred to as 
intermittency. The data indicate broader and broader tails for the distribution of 5u{l) at 
smaller and smaller scales, e.g., the kurtosis of the distribution, {6u{l)*) / {6u{iyy oc Z''^'^^''^, 
increases with decreasing I because ^4 < 2^2- The distribution of Su{l) is fatter for smaller 
I. The anomalous scaling is fundamentally caused by the fluctuations in the dissipation 
rate e^. Applying the refined similarity argument for homogeneous and isotropic turbulence 
(Kolmogorov 1962), the velocity difference over a separation I can be related to the dissipation 
rate e^, 

5u{l) ~ eJ/'Z^/^ (9) 
(note that the Kolmogorov's 1941 theory uses e.) The structure functions are then given by. 



{Sum « {ef)l'^'- 



(10) 
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Clearly, the departure from the linear scaling for the velocity difference comes from the 
statistics of the dissipation rate. Assuming (ef) oc (e.g., She and Leveque 1994), we have, 

Cp = p/3 + V3. (11) 

Developing a physical model for Tp that satisfies the experimental result for ({p) has been the 
main task of intermittency theories. Although we are mainly concerned with the distribution 
of ei, discussions of the structure functions are necessary because they are directly measurable 
in experiments and give important information and constraints on the ei distribution. We 
will use two intermittency models in our calculations: the log-normal model (Oboukhov 1962 
and Kolmogorov 1962) and the log- Poisson model by She and Leveque (1994). 



3.1. The log- normal model 

Oboukhov (1962) and Kolmogorov (1962) developed the first intermittency model. In 
this model, the distribution of ei is assumed to be log-normal (Kolmogorov 1962). A justifi- 
cation for this "natural" distribution for e; was given by Yaglom (1966). Imagine the cascade 
progress as successive eddy fragmentations from the integral scale L to the dissipation scale 
T). The statistics of the energy flux at an inertial scale I (or equivalently the dissipation rate 
e;) depends on the fragmentations before the scale is reached. The total number N of steps 
that lead to the scale / is proportional to iV ~ ln{L/l). Defining Xi = Ci/^-i as the ratio of 
the energy transfer rates at two successive fragmentation steps, the energy flux at the scale 
/ can be expressed in the ratios (see e.g., Monin & Yaglom 1975), 

ei = eLXiX2-XN (12) 

where is the dissipation rate at the integral scale (or the transfer flux at the largest scale). 
El — e. Due to the randomness in the fragmentation process, x^'s are stochastic variables. 
Assuming a self- similar fragmentation process, the distributions of Xj's are similar and 

N 

ln{ei/eL) = ^ in{xi) is expected to be Gaussian from the central limit theorem, 

i=l 

p{ei)dei = —==exp{- ^ ^ 2 )dHei/e) (13) 

V27r(7f 2(7; 

where the variance af is proportional to the number of steps af = iiln{L/l) with n being a 
parameter to be determined by experimental data and the (jf/2 term in the numerator in 
the exponential is to guarantee the mean (e^) is equal to the overall average dissipation rate 
e. This distribution will be used later to calculate the probability (eq. 4) for a region of a 
given size being in the distributed regime. 
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The scaling behavior of ei can be derived by integrating equation (13), 



(14) 



which gives Tp — —\np{p — 1). Prom eq (11), we have 

Cp = P/3 - :j^/^p(p - 3) 



(15) 



Therefore = 2 — (^g, which can be obtained from the results of experiments and simulations. 
It has been found that ji ~ 0.2 (Frisch 1995, Biskamp 2000). The relation (15) agrees with 
experiments quite well at small p but starts to exhibit deviation at p ^ 10 and gives an 
unrealistic maximum and turnover at p > 16, violating the requirement that the C(p)" P 
curve must be monotonic and concave (Frisch 1995). Simulations by Wang et al. (1996) 
suggest that this disagreement corresponds to the departure of the distribution for /n(e//e) 
from normal at scales close to the dissipation scale. They find that, at these scales, the 
distribution of ln{ei/e) shows a negative skewness, meaning that the log- normal distribution 
overestimates the probability in the very high e/ tail. However the distribution of agrees 
with log-normal very well in the inertial range away from the dissipation scale and the 
agreement is better and better for larger and larger scales (Wang et al. 1996). Fortunately 
the critical scale we are concerned with is well within the inertial range (see §4) and according 
to Fig. 6 in Wang et al. (1996) the log normal fit is very good at least up to the A — a 
tail. They also show that the fit gets better as the Reynolds number increases. More recent 
simulations by Yeung et al. (2006) with resolutions up to 2048"^ obtained similar results. The 
log normal distribution gives a very good fit to 4 — cr and only deviates by a factor of 2 at 
the 5 — (7 tail. They also find that the negative skewness gets closer to zero with increasing 
Reynolds number. The Reynolds number in SNe la is Re ~ 10^^ for typical velocity scale 10^ 
cm/s, length scale 10^ cm and viscosity 1 cm^/s. This is much larger than in all the current 
simulations. Therefore one may expect that the log-normal distribution probably applies 
even further out on the tail for the inertial scales of the turbulence in SNe la. However, the 
departure of the predicted ({p) — p curve from the experiments (with high Re) at p 10 
suggests that, even at huge Reynolds number, the log normal fit eventually breaks down at 
some large e; in the tail. Therefore we need to be careful when using the log-normal model. 
We will give more discussion on this point in the calculations given in §4. 

Another issue is that the distribution of ei has a physical cutoff in a realistic system due 
to the finite viscosity. Since the intermittency is stronger at smaller scales, the cutoff in the 
distribution of e; is probably larger for smaller I and obtains a maximum at the dissipation 
scale, 77. For Kolmogorov scaling, the cutoff in the distribution of is is given by ei?e^/^. 
Since Re ~ 10^^ in SNe la, this maximum dissipation rate is far beyond that required to 
break flames at density ^10^ g/cm^. Therefore ignoring this maximum cutoff does not 
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affect our result. However, the largest available dissipation rate at an inertial scale I is 
probably smaller than the cutoff in the distribution of and thus may affect the calculation 
for the cumulative probability defined in eq (4) if the cutoff in the distribution P{eiJ is close 
to or even smaller than K^^f- Since the log-normal model does not address the cutoff in the 
distribution of ei, we will neglect this potential effect in this model. 

On the other hand, the log-Poisson model we consider in the next section gives a maxi- 
mum dissipation rate at each inertial scale, corresponding to the strongest dissipativc struc- 
tures at that scale. In that model, a nonzero cumulative probability in eq (5) requires the 
lower limit K^ef in eq (4) be smaller than the maximum. 

3.2. The log-Poisson model 

A major success in the intcrmittency theory is the model by She and Lcvcquc (1994). In 
this model. She and Leveque studied the hierarchy of dissipation intensity in structures of size 
/ and, by invoking an unknown "hidden symmetry," they related the characteristic dissipation 
rates in structures of different intensity levels to the strongest dissipative structures. This 
relation gives a prediction of function of p, which only depends on the proprieties 

of the most intermittent structures. Assuming that the dissipation rate in regions of size I 
containing the most intense structures exhibits a scaling oc (see explanation in Appendix 
A) and the most intermittent structures are filamentary, corresponding to a codimension 
of 2, She and Leveque obtained a C(p) — P relation, which is in excellent agreement with 
experimental data. The "hidden symmetry" has been immediately interpreted as a log- 
Poisson process (Dubrulle 1994, She and Waymire 1995) in a multiplicative cascade model. 
In this section, we adopt the log-Poisson version of the She-Leveque model. The original 
presentation by She and Leveque (1994) is given in Appendix A. 

In a multiplicative model, the dissipative rates at two scales I2 and li {li > I2) are related 
by a multiplicative factor Wi-^i^, 

ei, = Wi,i,ei, (16) 

The average (Wi-^i^) is equal to unity since (e/^) = (e/j) = e. She and Waymire (1995) 
speculated that Wi^i-^ consists of two events. First is the amplification of the dissipation rate 
in the cascade, which tends to produce singular structures with ei^ oc {li/h)"' approaching 
infinity as I2 goes to 0. The meaning of 7 is discussed below. To ensure {Wi^i^) = 1, a second 
event is required to reduce Wi^i^. She and Waymire (1995) called this event the modulation- 
defects since it modulates the singular structures. The defects were assumed to be a discrete 
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Poisson process. Each of the defects decreases Wi^i^ by a factor of /3, thus 

W,j, = (17) 

if there are n defect events in the cascade. The number n of the events that occur in the 
cascade from the scale h to I2 obeys a Poisson distribution, 

P{n) = exp{-Xi,i,)^ (18) 
n! 

where Xi-^i^ is the mean number of the defect events in the cascade, which is expected to 

be proportional to the total number of the cascade steps, i.e., Aj^jj oc ln{li/l2). In fact, 

Xi^i^ can obtained by taking the average of eq. (17) and requiring {Wi^i^) = 1. Using the 

00 

identity ^ = exp{a), we get = exp{{l3 — l)A/^;2) ior the Poisson distribution eq 
(18), therefore, 

A,., = 2M1M (19) 

In this model, there is a largest dissipation rate at any scale. Clearly the largest dis- 
sipation rate is achieved if there is no defect, i.e., n = 0, in a cascade from the integral 
scale L to the scale / of interest, thus the largest dissipation rate is equal to ei^L/iy. This 
largest dissipation rate corresponds to e\°°^ in Appendix A. Similarly n = 1 gives the second 
strongest dissipative rate at a given scale, and so on. 

From eqs (16) and (17), we have, 

ln{ei^ /e) = ln{ei^ /e) + -flnih/k) + nln{P) (20) 

thus, using the Poisson distribution for n, the distribution for the dissipation rate at I2 can 
be derived from that at any scale h larger than I2. In particular, we consider deriving the 
distribution of ei at any scale I from the integral scale L. The distribution function of at 
the integral scale depends on how the energy is injected in the flow, thus is not universal and 
may vary from flow to flow. Therefore the function form cannot be specifled. However, there 
is a strong constraint for its width. Since ~ e, the distribution of /n(e/^/e) is expected to 
be very narrow around ln{tL/t) ~ and hence to be approximately a delta function. We 
denote the distribution of ln{eL/e) as PL{ln{eL/e)). It then follows from eqs (18) and (20) 
that 

P(ei)dei = V exp(-X)—PL(ln(ei/e) - -flniLll) - nlni(3))dlniei/e) (21) 

n=0 

where A = A^; = ^ln{L/l)/{l — P). Each term in eq (21) represents the contribution 
from dissipation structures of different levels, e.g., the n = term corresponds to the most 
intensive structures of size I. 



-16- 



To compare the model with experiments and obtain the parameters, we calculate the 
moments (e;) from the distribution eq. (21), 

oo 

(^f) '^Jexp{—X)^exp{px)PL{x — 'yln{L/l)—nln{P))dx 

n=0 

OO 

= {l/L)-^^exp{-X) E ^ / exp{x')P,{x')dx' (22) 
= Bpexp{-X{l - I3P)){1/L)-^P 

where we used a variable change x' — x—jln{L/l) —nln{f3) in the second step and the identity 

00 

^ — exp{a) in the third step. The coefficients Bp — J exp{px')PL{x')dx'; Bq — Bi — 1 
from the normalization of Pl and the requirement that (e^) — e respectively. 

The result eq (22) gives Tp = —jp + 7(1 — P'^)/ (1 — P), which is the same as (A8) in 
Appendix A, meaning that the "hidden symmetry" described in the appendix is equivalent 
to a log-Poisson process. The parameters 7 and P introduced here are identical to those 
described in the appendix, and thus have the physical meanings explained there, i.e., 7 can 
be interpreted as the exponent of the dissipation rate scaling in regions containing the most 
intermittent structures and P is related to the codimension C of the strongest dissipation 
structures, 7/(1 — /9) = C (see Appendix A for details). As discussed in the appendix, She 
and Levequc argued that 7 = 2/3 and P = 2/3 for C = 2 corresponding to filamentary 
dissipation structures in incompressible turbulence. This results in (p as a function of p 
that agrees with the experiments with an accuracy of 1%, implying eq (21) provides a good 
distribution for e^. The She-Leveque formulation has been extended to supersonic turbulence 
(Boldyrev et al. 2002) and MHD turbulence (MuUer & Biskamp 2000) where the dissipation 
structures are dissipation sheets and the current sheets, respectively. For these 2 dimensional 
dissipation structures, the codimension C = 1 and P ~ 1/3. In next section we use the log- 
Poisson distribution (eq 21) in our calculations for the cumulative probability in eq 4. We 
will take 7 = 2/3 and consider both filaments {P = 2/3) and sheets (/3 = 1/3) as the most 
intermittent dissipation structures. 
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4. Results 

4.1. The log- normal model 

We are ready to calculate the probability -P(e/^ > K^ef) using the distributions P{ei) 
given in §3. The calculation is straightforward for the log-normal distribution eq. (13), 



ln{K'iti/l) 



oo 



= / -^^expi^-x"^ /2)dx (23) 

- 2^^J^K V2ai^ + 2V2) 

where e ~ U^/L, af^ = nln{L/lc) ~ 0.2/n(L//c) and erfc{x) is the complementary error 
function. Using e/ and Ic given in Table 1, we calculated the probability as a function of the 
density assuming different values for the characteristic velocity (U) and length (L) scales. 
For example, if [/ = 100 km/s and L = 100 km, the probabilities are 0.5 x 10"^°, 0.042, 0.9 
and 1 at p = 10^, 5 x 10^, 3 x 10^ and 10^ g/cm^ respectively if = 1. It is interesting to note 
that, at a density 3 x 10^ g/cm^, 10% of the local regions of the critical size are still in the 
flamelet regime, although generally the flame has reached the distributed regime according 
to the mean criterion eq (2). With i?j ~ L = 100 km, at the 4 densities above from high to 
low, the corresponding numbers Nic = AirR'j/l^ of regions of the critical size that cover the 
flame front are 3 x 10^°, 7 x 10^ 5 x 10^ and 3 x lOl Multiplying P(q^ > K^tf) with A^;^ (eq. 
5), we see that there is already one region of size Ic in the distributed regime at a density 
10^ g/cm^. Recalling that, according to the criterion eq (2), the DDT does not occur until 
the density decreases to 4 x 10^ g/cm^, we find that in this case the intermittency effect may 
increase the transition density by more than a factor of 2. 

We point out the cumulative probabihty calculated from eq (23) at density 10^ g/cm^ 
in the example above comes from a little beyond the 6 — 0" tail of the distribution for ln{ei/e) 
and we need to check whether the log-normal distribution there is a good approximation. 
As discussed in §3.1, numerical simulations have shown that for a scale / in the inertial 
range, the distribution of e; is well approximated by log-normal up to the 5 — o" tail (Yeung 
et al. 2006). Assuming the Kolmogorov scaling, the dissipation scale in Type la SNe is 
7] = LRe^^/^ ~ 10~^ cm for L ~ 10^ cm and the Reynolds number Re ~ 10^^. The critical 
scale Zc ~ 10^ — 10^ cm of interest here is well between the integral scale and the dissipation 
scale, thus we expect that the distribution for ei^ is close to log-normal at least up to the 
~ 5 — (7 tail. The question is then whether the good fit extends further. Wang et al. 
(1996) and Yeung et al. (2006) found that the log- normal fit is better for larger Reynolds 
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number, thus it is expected that the log-normal approximation probably applies to higher 
on the tail than 5 — 0". As argued in §3.1, the log normal approximation eventually fails 
somewhere in the extreme tail even at high Reynolds number. To know exactly how far the 
log-normal fit extends, numerical simulations with much higher resolution are needed. We 
have to be careful about the validity of the log-normal approximation in the far tail because 
it overestimates the probabihty distribution for ei once it breaks down and in that case eq 
(23) overestimates P{ei^ > K^^f)- 

Due to the complication of the validity of the log-normal distribution at the far tail, we 
consider two extreme cases and give the upper and lower limits for the transition density. 
First, we ignore the departure from log-normal and evaluate the density at which Ni^xP[ei^ > 
K^ef) = 1 using eq (23) for P{ei^ > K^ef) and Table 1 for e/ and Ic with different parameters 
U and L. We will denote this density as pln with the subscript LN standing for log-normal. 
Interpolation was used to obtain e/ and Ic not tabulated in Table 1. If the distribution of 
6; is exactly log-normal as given by eq (13), then p£,7v is the predicted transition density 
for the DDT with the intcrmittcncy taken into account. On the other hand, if the log- 
normal distribution overestimates the probability at the high tail, eq (23) overestimates 
the cumulative probability and p^^ is the upper limit for ptr- We give p^jv for different 
parameters U and L in the second line of Table 2. 

In the other extreme, we assume that the log-normal distribution fails to fit the distri- 
bution of ei beyond the 5 — a tail. This gives a lower limit for the transition density since 
numerical simulations have shown that the log-normal fit applies at least to 5 — o". In this 
case, we keep track of the integral limit in the second line of eq (23) at pln, which tells us 
which part of the tail of the distribution gives the main contribution to P{ei^ > K^ef) at that 
density. If the integral limit is smaller than 5, the contribution to the cumulative probability 
is from within 5 — a and vice versa. We calculate the density at which the integral limit is 
equal to 5 and denote this density as p^^- Since the integral limit is a decreasing function of 
the density, if pln < P5a, the contribution to P{ei^ > K^^f) at density pln is from within 
5 — 0". In this case, the cumulative probability calculated from eq (23) is valid and pln is 
a good estimate for the transition density. Otherwise if p^jy > P5a, the contribution to the 
probability is from beyond the 5 — 0" tail, eq (23) overestimates it and thus p^N overestimates 
the transition density ptr- In this case, p^o- gives a lower limit for the transition density be- 
cause at P50-, we have Ni^ x P{€i^ > K^ef) 3> 1 using eq. (23) which apphes for p < ps^-. 
Therefore if the log-normal fit fails just beyond 5 — o", we have a lower limit for the transition 
density, min{pLN, p^a)- We give this lower limit in the 3rd line of Table 2. 

Similar calculations can be done for the K — % case. The results of the upper and 
lower limits for the transition density in the K — 9) case are given in parenthesis in Table 
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2. Comparing with predictions from the mean criterion eq (2) (the first hne in Table 2), the 
log-normal model predicts that the intermittency effect increases the transition density by a 
factor of 2-3 for all the cases we hst in Table 2. 

We evaluate the probabihty of the existence of a region of size 10^ cm with a rms velocity 
10^ cm/s, required for the DDT by Lisewski et al. (2000), and compare with the numerical 
results of Ropke (2007). Using the log-normal distribution for e/, we find that the requirement 
-^(^106 cm > 10^^ cm^/s^) requires conditions from the extreme tail of the distribution. The 
likelihood is completely negligible (~ 10^^°) if the velocity U at the integral length scale is 
less than ~ 10^ cm/s. Only if U is larger than 5 x 10^ cm/s is the probability appreciably 
larger so that the required region might be available. For example, if U = 5 x 10^ cm/s at 
L — 10^ cm, the probability is ~ 10~^^. This is still too small to guarantee the existence of 
a region as required by Lisewski et al. (2000). The number of available candidate regions of 
size 10^ cm around the flame front is probably smaller than 10^ — 10^, assuming the flame 
front radius is ~ 10^ — 10^ cm. This result agrees with the conclusion of Ropke (2007) that 
the existence of a region as required by Lisewski et al. (2000) is rare. To ensure such a 
region, the velocity at the integral scale has to be larger than 10^ cm/s, which is probably 
impossible as discussed in §2. 

We also carry out a calculation for ptr based on the criterion of Woosley (2007) taking 
into account the effect of intermittency. Using the log-normal distribution to calculate the 
cumulative probability in eq. (7), we find that no regions of critical size that meet Woosley's 
criterion appear at density above 10^ g/cm^. We cannot give an exact predicted transition 
density for this model because we do not have data at densities below 10'' g/cm^ for relevant 
quantities hsted in Table 1. Note that Woosley (2007) obtained a transition density around 
10^ g/cm^ under the assumption that a region of size 10^ cm with rms velocity of 10^ cm/s is 
available. Prom our estimate above and the result in Ropke (2007), the probability that such 
a region exists is small, therefore it is appropriate to take the transition density predicted in 
Woosley (2007) as an upper limit for his DDT criterion. 
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4.2. The log-Poisson model 

We next consider the log-Poisson model. Using the distribution eq (21), we have, 

oo oo 

P{€i^ > Khf) = exp{-K) E f / Pl{x - iln{L/Q - nln{(5))dx 

n=0 ln{K^ef/t) 
oo oo 

= exp(-Ae)Ef / PL{x)dx (24) 

n=0 ln{K^tf/t)--iln{L/lc)-nln{li) 

oo 

= exp{-Xc) E 

n=0 

where Ac = ^ln{L/lc)/{l — (3) and the integrals in the second line are denoted as F„ for 
convenience. Note that the integral lower limit ln{K^ef/e) — ^ln{L/lc) — nln{P) increases 
with n because P < I, therefore F„ is a decreasing function of n. An exact calculation for the 
cumulative probability is impossible because of the unspecified function Pl. We will neglect 
all the n > 1 terms and only keep the n = term in our calculation, i.e., we only include 
the contribution of the most intensive structures at scale l^. Obviously, this approximation 
gives a lower limit for the probability and the transition density we obtain will also be a 
lower limit. We will show that the criterion for DDT obtained from this approximation is 
exact if Pl is a delta function. 

The contribution from n = is exp{—Xc.) = Fq. For 7 = 2/3 and P = 2/3, it 

is equal to (h/LfFo and Ni^ x P(e,, > Khf) > A-K{Rf /l^f{lj Lf Fo = A'ni.Rf/LfFo. Since 
the size of the flame region Rf > L, it means that the number of regions which are larger 
than the critical size and in the distributed regime is ~ AttFq. Since the distribution Pl{x) 
is probably strongly concentrated at a; = 0, the sufficient and almost necessary condition for 
Fq ~ 1 is that the integral limit ln{K^€f/e) — ^ln{L/lc) < 0, or equivalently, 

e- > {IjLf/^Khf (25) 

which is a convenient criterion for the DDT in the log-Poisson model. Note this criterion 
is much weaker than the mean criterion eq (2). Once the condition is satisfied, at least one 
region of critical size that covers the fiame enters the distributed regime due to the most 
intense stretching strength available at scale Ic- 

As mentioned in §3.2, if the dissipation structures are 2-dimensional, (5 — 1/3. In 
that case, the contribution from the n = term is {lc/L)Fo and Ni^ x P(e;^ > K^ef) > 
'i7r{Rf/lc){Rf/L)Fo, which is much larger than 1 if Fq 1. Therefore, the criterion eq (25) 
is a sufficient condition for the case with sheet-like dissipation structures such as in MHD 
turbulence or highly compressible turbulence. 

We have neglected the n > 1 terms in eq. (24), the contribution of which depends on 
how rapidly Pl{x) decreases with x > 0. We consider the extreme example where Pl is a 
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delta function. In this case, before the condition eq (25) is met, F„ = for any n thus the 
cumulative probability is zero. When the condition is just satisfied as the density decreases, 
only the n = term contributes and all the n > 1 terms are still zero, i.e., the most intensive 
(n = 0) structures at Ic can stretch a local flame into the distributed regime while all the less 
intensive structures (n > 1) still cannot. From the calculation above, we see that in this case 
once the n = term contributes, at least one region around the flame front experiences the 
largest stretching rate and enters the distributed regime. Therefore, if Pl is a delta function, 
eq (25) is both the necessary and the sufficient condition. This is true for both P — 2/3 and 
P — 1/3. If Pl is not a delta function, the tail of Pl gives rise to the possibihty that the 
distributed regime can emerge in a local region of critical size before the condition eq (25) 
is met. This could lead to an even weaker condition than eq (25). Since we expect that cl 
can only vary within a factor of a few, the condition can be weaker only by a factor of a few. 
Because the r.h.s of the condition (23), especially e/, depends on the density very sensitively, 
this would not increase the predicted ptr considerably. 

The condition eq (25) can be easily applied to calculate the transition density using e/ 
and Ic in Table 1. For example, we get Ptr = 8.7 X 10^ g/cm^ for U = 100 km/s and L = 100 
km if K = 1. This result is consistent with that from the log- normal model and is also about 
a factor of 2 larger than the prediction by the mean criterion eq (2). 

The transition density predicted by the log-Poisson model with different parameters for 
turbulence is given in the 4th line of Table 2. Again the numbers in parenthesis are for 
K = 8. The results arc consistent with those from the log- normal model and are at least 
2 — 3 times larger than from the criterion eq (2). 

Again we consider the possibility that there exists a region of size 10^ cm with rms 
turbulent velocity of 10^ cm/s required for DDT by Lisewski et al. (2000). The probability 
P{^io^ cm > 10^^ cm^/s^) depends on Pl, the probability distribution of the dissipation rate 
at the integral scale L. Since Pl is probably not universal and is flow- dependent, the log- 
Poisson model cannot give an exact estimate for the probability. Here we assume Pl is a delta 
function and see under what condition it is possible to find a required region. We find that 
the necessary and sufficient condition to have such a region is that U > 10^^/^(L/cm)^/^ cm/s. 
For L ~ 10^ cm, U has to be larger than 10^ cm/s. This can be understood from the fact that, 
in the log-Poisson model, the available kinetic energy in the most intermittent structures for 
dissipation is assumed to be the kinetic energy at the integral scale (see Appendix A). Since 
U > 10^ cm/s is probably not achievable, it is rare that a region as required by Lisewski et 
al. (2000) exists, again in agreement with Ropke (2007). 

Using the same calculation that leads to eq (25), we obtain a criterion for the DDT 



-22- 



model by Woosley (2007) accounting for the intermittency effect, 

e>{ljLf/\lJljfef. (26) 

This is weaker than the corresponding mean criterion eq (6) by a factor of {Ic/Ly^^, meaning 
that intermittency increases the transition density. In comparison with eq (25), the condi- 
tion is stronger and thus gives a smaller transition density than that for the KOW and NW 
model with intermittency included. At smaller density, the critical length Ic is larger and 
the factor [Ic/LY^^, representing the intermittency effect, is closer to unity. This imphes 
that intermittency gives a weaker effect on the the transition density for the Woosley (2007) 
criterion than for that by KOW and NW. Using Table 1, we again find that the condition 
eq (26) is not satisfied at densities above 10'^ g/cm^ for the five cases listed in Table 2, i.e., 
the predicted transition density is still below 10^ g/cm^ after including intermittency (see 
discussion in §4.1). 



In summary, intermittency can considerably enhance the onset of the distributed fiame 
regime and hence increase the transition density in the DDT model of KOW and NW. Both 
the intermittency models we consider here predict a transition density 2 — 3 times larger than 
from the criterion using the mean dissipation rate. This factor of 2 — 3 brings the transition 
density to be in disagreement with the observational constraints for turbulent velocity larger 
than U = 10^ cm/s in the case K — 1. We discuss the imphcations of this result in the next 
section. We also find that existence of regions of size 10^ cm with velocity 10^ cm/s is rare, 
in agreement with the numerical result of Ropke (2007). The strong DDT criterion given by 
Woosley (2007) gives a transition density below ptr = 10^ g/cm^ even when intermittency is 
included. We expect that the intermittency effect is weaker for stronger DDT criteria. 

5. Conclusion and Discussion 

We have studied the effect of intermittency on the transition from the flamelet regime 
to the distributed regime in Type la SNe, and hence on the transition density for the DDT 
model by KOW and NW. In their model, the detonation occurs via the Zel'dovich mechanism 
that requires a nearly isothermal region larger than a critical size to drive a sufficiently strong 
supersonic shock. KOW and NW assumed that the almost isothermal mixture of fuel and 
ash can be produced once turbulence is strong enough to get the flame into the distributed 
regime. The DDT is assumed by KOW and NW to occur when the average flow gets into 
the distributed regime. We argue that the sufficient condition for the DDT is that there is 
one region that is larger than the critical size and in the distributed regime. 
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The intermittency in turbulence, as a result of the spatial inhomogeneity of the dis- 
sipation rate, gives rise to regions with strong local turbulent strength that can force the 
flame into the distributed regime earlier than elsewhere. Therefore the transition from the 
flamelct regime to the distributed regime is not spatially smooth, but intermittent. At early 
time when the density in the white dwarf is large, the flame has a large speed and a small 
width and thus resists being efficiently stretched and broken by the turbulence. At the same 
time, the critical size is very small. This has two effects that tend to make an early DDT 
likely. First, the intermittency of turbulence tells us that the probability of finding extremely 
strong stretching within a smaller critical size is larger. Second, there are more regions of 
smaller sizes available. Therefore it is possible that the DDT is triggered at a small "spot" 
when the density is larger than needed for the average flow to enter the distributed regime. 
As we pointed out in the Introduction, the critical size as a function of the density plays an 
important role in determining the transition density for the DDT in our calculations. 

We used two analytical intermittency models to statistically investigate when the flrst 
region appears which is both larger than the critical size and in the distributed region. This is 
assumed to be the time when the DDT occurs by KOW and NW. We found that, for various 
parameters for the intensity and length scales, DDT occurs at a transition density at least 
2 — 3 times larger than the density at which the average flow enters the distributed regime. 
The transition density has been determined empirically by invoking it as a free parameter 
in spherically-symmetric models and then computing models that best match the observed 
multicolor hght curve shapes and magnitudes (Hoflich & Khokhlov 1996). Recognizing that 
the spherical models are oversimplified, they do give some guidance to the empirical con- 
straints on the density at which DDT occurs. Hoflich (1995) used this procedure to flt 
observations of the Branch core normal SN 1994D and preferred a value of the transition 
density of 2 x 10^ gm cm^'^. Hoflich, Khokhlov & Wheeler (1995) explored a range of transi- 
tion densities in the context of pulsating delayed detonation models and favored densities in 
the range 0.8 — 2.2 x 10^ gm cm'^. Dominguez, Hoflich & Straniero (2001) adopted 2.3 x 10^ 
gm cm~^. Allowing for an uncertainty of a factor of 2, the predicted transition densities by 
the mean criterion are consistent with 2 x 10^ gm cm~^ as favored by the observations in all 
the cases except that with e = 10^^ cm^/s^ and K — 1 (Table 2). With the intermittency 
effects we have examined here, the transition density would be a factor of 2 — 3 higher. If 
i^T = 1, all the predicted ptr are larger than 2 x 10^ gm cm^^ by at least a factor of 2 except 
the case with U = 10^ cm/s. The predicted transition density with i^T = 8 is 2 — 3 times 
smaller than from K = 1. From Table 2, the predicted ptr for the intermittency models with 
K = 8 agree with the observations within a factor of 2 except the case with a large velocity 
scale U — 10^ cm/s at the integral length scale. To avoid discrepancy with the observations, 
our result indicates several possibilities. 
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1. The large scale motions caused by Rayleigh- Taylor instability freeze out due to the 
overall expansion of the star (Khokhlov 1995). The freezeout effect has to be efficient 
enough so that the developed part of the flow has a velocity scale of ^ 10^ cm/s (see 
Table 2). 

2. The flame is very robust. To break the flame, the local Gibson scale has to be at least 

— 512 times smaller than the flame thickness. In this case, the predicted transition 
density is 2 — 3 times smaller than from K — 1. 

3. There is not enough time for the buoyancy-driven turbulence to fully develop down 
to the critical size before the predicted density for the DDT by our intermittency 
models is attained, thus motions at scales below the critical size are either absent or 
non- intermittent . 

4. The DDT does not occur immediately after a region of the critical size enters the 
distributed regime. It may take some time for turbulence to help mix the region and 
make it nearly isothermal. However, the time scale for turbulence to mix a region of 
the critical size in the distributed regime is very small, ^ 10~^ s at densities larger 
than 3 x 10^ g/cm^. It is unlikely that the density drops much in such a short timescale. 

5. Having a large enough region entering the distributed regime is not a sufficient condi- 
tion for detonation. As mentioned in the Introduction, there are several uncertainties 
in the simple model by KOW and NW assuming flame quenching, entering the dis- 
tributed regime and the DDT all occur simultaneously. 

Our result alludes to the possibility that the criterion by KOW and NW is too weak for 
the DDT, supporting the claim of Lisewski et al (2000) and Woosley (2007) that just entering 
the distributed regime is not sufficient for the DDT. We have shown that their criteria for 
the DDT are much stronger than just entering the distributed regime. We also studied the 
intermittency effect on their conditions for the DDT. We find that the existence of a region 
of size 10^ cm with rms turbulent velocity of 10^ cm/s required by Lisewski et al (2007) for 
a DDT is rare, consistent with numerical results of Ropke (2007). We have also examined 
the intermittency effect on the transition density for the DDT criterion by Woosley (2007). 
We find that the effect is weaker for the stronger criterion and does not increase ptr to above 
10^ g/cm^. Woosley (2007) obtained ptr around 10'' g/cm^ because of the assumption of 
strong turbulence velocity of 10^ cm/s in a region of size 10^ cm. Since the existence of 
such a region is rare, Woosley (2007) may considerably overestimate the transition density. 
This may imply that the condition for DDT by Lisewski et al. (2000) and Woosley (2007) is 
too strong and predicts a transition density smaller than that empirically determined from 
observations. 
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We have only studied the DDT in white dwarfs with the initial chemical composition 
of half carbon and half oxygen. In a white dwarf with more carbon, the nuclear timescale 
is smaller thus a stronger turbulent intensity is needed to break the flame. This results in a 
smaller transition density using the mean criterion (NW). On the other hand, the critical size 
for detonation in such a white dwarf is smaller, therefore the intermittency effect could be 
more efficient in increasing the transition density. We were not able to perform a calculation 
for the chemical composition with the carbon abundance larger than 0.5 due to the lack of 
sufficient data for the critical length scale in this case. 

We point out that the intermittency models we used were originally proposed for ho- 
mogeneous and isotropic turbulence. Turbulence in SNe la is stratified and may not achieve 
homogeneity and isotropy at very small scales even if the turbulence is developed at these 
scales. The effect of the departure from homogeneity and isotropy on the predicted transition 
density is out of the scope of this paper. 

We thank Elaine Oran for useful discussions. This research was supported in part by 
NSF Grant AST-0707769 (LP, JCW) and by NASA Astrophysics Theory Program Grant 
NAG5-13280 (JS). 



A. The She-Leveque Model 

In their original paper. She and Leveque (1994) start with the moments of the distribu- 
tion P{ei) and use the ratios of two successive moments, ef^ — (ef"*"^) / (ef ) oc I'^p+^-'^p^ to char- 
acterize a hierarchy of dissipative structures. This ratio can be written as e^^ — J eiQp{ei)dei 
where Qp{ei) = e^P{ei)/ J e^P{ei)dei. For a typical distribution P{ei) that decreases mono- 
tonically and faster than any power law at large e;, Qp{ei) strongly peaks around e[^^ for 
large p. Clearly el^'* increases with p, and e;"^"* = lim„^oo(ef'''^)/(ef) corresponds to the most 
intense dissipative structures at scale /. These strongest dissipative structures are the origin 
of the anomalous scaling and the scaling of e[°°^ with /, 

eS°°) ~ (Al) 

is of fundamental importance. To determine the parameter 7, we can dimensionally write 
e\°°^ as an energy scale divided by a time scale ti. She and Leveque argued that for the 
most intermittent structures this energy scale is the largest available kinetic energy (which 
is ~ "^^rms! independent of /) and assumed that ti exhibits a regular Kolmogorov scaling 
ti ~ therefore 7 = 2/3. From eq (Al), we have r^+i = — 7 for p ^ 00 or 



(A2) 
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where the constant C has a physical interpretation as the codimension of the most intermit- 
tent structures. Eq (Al) means that the dissipation rate in a region of size / that encloses the 
most intensive structures scales with / as e; oc . When calculating (ef ) at p — > cxo, we need 
to consider the possibility of a point finding itself within a distance / to the most intermittent 
structures, which is proportional to /(^""^^ where D is the dimension of the system and d is 
the dimension of the most intensive structures (Frisch 1995). As p — > cxo, the contribution to 
(ef ) is dominated by the most intermittent structures, therefore (ef) oc /-TP+^-fi, Compar- 
ing with eq (A2), we find that C corresponds to the codimension of the most intermittent 
structures, C — D — d. 

In order to determine the entire hierarchy of the dissipative structures. She and Lcvcque 
argued that the intensity e^^"^^^ of the dissipative structures at level p -\- 1 depends only on 
their immediate precursor, the structures of level p, from which the level p -\- 1 structures 
directly develop, and on the most intensive structures, where the structures of all orders 
tend to end up. Based on this argument, they made an assumption about the hierarchy of 
the dissipation rates, 

e'r'-M^fYi^'r'f-' (A3) 

where the coefficients Ap are independent of / but may be flow-dependent and non-universal. 
The parameter /3 will be completely fixed by 7 and the codimension C. According to She and 
Leveque, this relation corresponds to a mysterious symmetry of the Navier-Stokes equation, 
termed "the hidden symmetry". 

To derive Tp from equation (A3), it is convenient to define a new variable, 

vr/ = ei/ef°^ (A4) 



which was introduced by DubruUe (1994). Clearly e^/ef""^ = (7rr^>/«> ^om the definition 

assumed by She 



of ef* . Then the "hidden symmetry" assumed by She and Leveque becomes. 



This recursion relation is solved by. 



(A5) 



^"^^'^ = Cp{7:ir (A6) 



p 

n 

n=0 



where Co = 1 and Cp = fl ^n" for p > 0. Eq. (A6) gives, 

= (A7) 
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where Bq = Bi = 1 and Bp = Yl Cn for p > 1. Noting that (tt;) = e/ef^' oc (e^ and 

71=1 

(7rf) = (ef)(e(°°))-f, eq (A7) gives (e[) oc (e;~))f-(i-/3^)/(i-/3). Using eq (Al), we have 

fl - 0^) 

Tp = -IP + T^Y^y^ ^^S) 

The parameter j3 is determined by the asymptotic behavior of Tp at the p — > cxd, eq (A2). 
Letting p — > oo in eq (A8) and comparing with eq (A2), we find that 7/(1-/3) = C — D — d. 
Since the most intermittent structures in 3D incompressible turbulence are filamentary, we 
have d — 1 and C — 2, thus P — 1 — 7/2 = 2/3 for 7 = 2/3. Finally we arrive at the 
celebrated She-Leveque formulae, 

T^ = -2p/3 + 2(l-(2/3f) (A9) 

and 

Cp = p/9 + 2(l-(2/3f/3) (AlO) 

which agrees with the experimental result with an accuracy of 1%. Note that this result 
is consistent with the Kolmogorov's exact result for the third order structure function, i.e., 
^3 = 1. If the most intense structures are two dimensional, e.g., the dissipation sheets in 
compressible flows (Boldyrev et al. 2002) or the current sheets in MHD (MuUer & Biskamp 
2000), /?= 1/3. 
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Table 1: The laminar flame speed, the flame thickness and the critical length for a white 
dwarf with half carbon and half oxygen 



p (10^ g/cm3) 


Si (10^ cm/s) 


If (cm) 


ef (IQis cmVs3) 


/c(cm) 


2 


75.8 


9.35(-5) 


4.66(9) 


7(1) 


0.5 


18.1 


9.46(-4) 


6.27(6) 




0.1 


2.33 


2.75 (-2) 


4.60(2) 


2(2) 


0.05 


0.599 


5.19 (-1) 


0.414 


1.3(3)^ 


0.03 


0.26^ 


1.78^ 


0.98 (-2) 


5(3) 


0.01 


4.72(-2) 


4.22 


2.59(-5) 


2(5) 



''Read from Fig. 6 in KOW. 
''Read from Fig. 7 in KOW. 

Note. — The values of Si and If are mainly taken from Table 3 of Timmes and Woosley 1992. The value 

of Ic is mainly taken from NW. We also include numbers (marked) from KOW because their results are very 
similar to NW despite the difference of details in the two models. Numbers in parentheses are powers of 10. 
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Table 2: The ])re(licto(l transition densities /-;/,. (in nnit of lO" g/cnr^) for various models 



Cases * 


A 


B 


C 


D 


E 


mean criterion 


6.9(3.3) 


5.5(2.3) 


4.1(1.5) 


3.0(1.0) 


2.0 (-) 


log- normal ^ 


27(10) 


11(4.9) 


10(4.4) 


9.4 (4.0) 


4.3 (-) 


log- normal ^ 


23(7.6) 


9.7(4.9) 


8.4(3.9) 


7.1 (3.0) 


4.3 (H 


log- Poisson 


24(8.6) 


9.5(4.3) 


8.7(3.8) 


7.9(3.3) 


3.8 (-) 



"A: U = 10^ cm/s, L = 10^ cm, e = IQi'^cmVs^; 
B:U= 10^ cm/s, L = 10^ cm, e = lO^^cmVs^ 
C:U = 10^ cm/s, L = 10^ cm, e = lO^cm^/s^; 
D:U= 10^ cm/s, L = 10* cm, e = lO^-'^cmVs^; 
E:U = 10^ cm/s, L = lO" cm, e = lO^^cm^/s-'' 

''The predicted transition density assuming a perfect fit of the distribution P{eiJ by log- normal. This is the 
upper limit for ptr since the log-normal approximation may break down and overestimate the distribution 
at the far tail. 

'^The lower limit for the transition density assuming the log-normal approximation applies only up to 5 — a. 

Note. — The numbers in parentheses are the results predicted if the Gibson scale has to be 512 times 

smaller than the flame thickness for the transition to the distributed regime. In case E, ( — ) indicates that 
the transition density is smaller than 10^ g/cm^, which cannot be well estimated since we only have data to 
10'^ g/cm^ in Table 1. 



